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Abstract 

We construct a class of matrix models, where super symmetry (SUSY) is spontaneously broken at 
the matrix size N infinite. The models are obtained by dimensional reduction of matrix-valued 
SUSY quantum mechanics. The potential of the models is slowly varying, and the large- N limit 
is taken with the slowly varying limit. 

First, we explain our formalism, introducing an external field to detect spontaneous SUSY 
breaking, analogously to ordinary (bosonic) symmetry breaking. It is observed that SUSY is 
possibly broken even in systems in less than one-dimension, for example, discretized quantum 
mechanics with a finite number of discretized time steps. Then, we consider spontaneous SUSY 
breaking in the SUSY matrix models with slowly varying potential, where the external field is 
turned off after the large-iV and slowly varying limit, analogously to the thermodynamic limit 
in statistical systems. On the other hand, without taking the slowly varying limit, in the SUSY 
matrix model with a double-well potential whose SUSY is broken due to instantons for finite 
N, a number of supersymmetric behavior is explicitly seen at large N. It convinces us that the 
instanton effect disappears and the SUSY gets restored in the large- N limit. 



1 Introduction 



Motivated by possibility of supersymmetry (SUSY) breaking in large-iV matrix-model 
formulations of superstring theory [HE], [3] (in particular [2]), we presented some concrete 
models, where SUSY is preserved in any finite size N of the matrix variables, but gets 
spontaneously broken at infinite N [I]. They are in the form of matrix models coupled 
to some supersymmetric field theories. The supersymmetric field theories alone undergo 
spontaneous SUSY breaking, however the couplings to the matrix model sector supply 
some effect to prevent the breaking. Thus, the models for finite N preserve the SUSY. As 
N gets large, the effect diminishes, and eventually nothing prevents the SUSY breaking 
in the field theory sector at N infinite. In this mechanism, the models exhibit the SUSY 
breaking only at the large- N limit. 

From our motivation, it is certainly desirable to construct models possessing the above 
property within the framework of a matrix model, not relying on SUSY breaking in 
the field theory sector. In this paper, we construct a class of such desirable large-iV 
SUSY matrix models. They are obtained by dimensional reduction of matrix-valued 
supersymmetric quantum mechanics with Af = 2 SUSY. It is expected that they play 
some role to reveal the essence of SUSY breaking of nonperturbative formulation of string 
theory via matrix models. 

In the next section, analogously to the situation of ordinary spontaneous symmetry 
breaking, we introduce an external field to choose one of degenerate broken vacua to 
detect spontaneous SUSY breaking in SUSY quantum mechanics. The external field 
plays the same role as the magnetic field in the Ising model introduced to detect the 
magnetization. To introduce a suitable external field for the supersymmetric system, 
we deform the boundary condition for fermions from the periodic boundary condition 
(PBC) to a twisted boundary condition (TBC) with twist a, which can be regarded as 
the external field. If a supersymmetric system undergoes spontaneous SUSY breaking, 
the partition function with the PBC for all the fields, Z PBC , which usually corresponds to 
the Witten index [5], is expected to vanish. Then, the expectation values of observables, 
which are normalized by Zpbc, would be ill-defined or indefinite. By introducing the twist, 
the partition function is regularized and the expectation values become well-defined. It 
is an interesting aspect of our external field for SUSY breaking, which is not seen in 
spontaneous breaking of ordinary (bosonic) symmetry. 

Differently from quantum field theory whose degrees of freedom are infinite, one- 
dimensional SUSY quantum mechanics has finite degrees of freedom, and thus the super- 
selection rule does not hold in general. Here, introducing the twist a enables us to show 
that, concerning the expectation values of the auxiliary fields, which are order parameters 
of the SUSY breaking, the superselection rule works and the expectation values at the 
one broken ground state are well-defined. On the other hand, we see that the expectation 
value of the scalar field becomes ill-defined as the external field a is turned off. 

Furthermore, it is straightforward to apply our argument to systems in less than one- 
dimension, for example discretized SUSY quantum mechanics with a finite number of 
discretized time steps. Spontaneous SUSY breaking is observed even in such simple sys- 
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terns with lower degrees of freedom. Also, we give some argument that an analog of the 
Mermin-Wagner-Coleman theorem [HI Ej does not hold for SUSY. Thus, cooperative phe- 
nomena are not essential to cause spontaneous SUSY breaking, which makes a difference 
from spontaneous breaking of the ordinary (bosonic) symmetry 0. As for the superselec- 
tion rule, the situation is similar to the SUSY quantum mechanics. However, for large- N 
SUSY matrix models discussed from section [3J we can see that the superselection rule 
does hold in general thanks to the large- N limit. 

In section [3J we define a matrix model analog of the discretized SUSY quantum me- 
chanics with general superpotential, and give an overview of its general aspects in the 
large- N limit. We mainly discuss the simplest case with the number of the time step 
T = 1. To realize SUSY breaking at large N, we put a parameter e into the super- 
potential. 1/e represents the length scale of the scalar in which the derivative of the 
superpotential W(</>) varies by the amount of 0(1). Thus, as e gets smaller, the potential 
becomes more slowly varying. If we take the large- N limit with the slowly varying limit 
e — > 0, the system is shown to exhibit spontaneous SUSY breaking, a is an infinitesimal 
external field which slightly breaks the SUSY, and we observe whether its breaking ef- 
fect remains after the limit a — * following the large- N and slowly varying limit, which 
is analogous to detect spontaneous symmetry breaking in the thermodynamic limit of 
a statistical system. As a deformation of the model, when the twist a is changed to a 
field-dependent function, the SUSY is also spontaneously broken in the same limit. The 
eigenvalue distribution of does not qualitatively change by the deformation, but its 
precise form is modified. Since the deformation effect remains even after turning off the 
external field, we see that the system in the large- N limit with e — > is somewhat sensitive 
to the deformation of the external field. Thus, as a complete definition of the model, we 
should specify not only the matrix-model action but also the form of the external field. 
This is reasonable because the twist modifies the boundary condition of the fermions, 
which is a part of the definition of the model even if the twist is turned off eventually. 

In section HJ we discuss one example of the models discussed in section [3J where the 
scalar potential is in a double-well type (W'((fi) is quadratic). Ref. [5] discusses the two- 
dimensional Wess-Zumino model with such a double-well potential. For finite spatial 
volume, the instanton effect causes SUSY breaking, while the effect ceases in the infinite 
volume limit to restore the SUSY. Analogously, in our matrix model case, SUSY will 
be spontaneously broken due to some sort of instanton effect for finite N. However, if 
we take the large- N limit with e kept finite (not slowly varying), the instanton effect is 
expected to disappear and the SUSY will be restored. By examining large-iV solutions 
of the matrix model, we observe a number of supersymmetric behavior that strongly 
indicates the restoration of the SUSY at the large- N limit. To best of our knowledge, this 
is the first example of a matrix model with the restoration of SUSY at large N, although 
Ref. discusses it for a large- N vector model. 

In the SUSY matrix model with a double-well potential, we find an intriguing fact 
that the free energy does not depend on the ratio of the numbers of the eigenvalues filled 

1 It has already been pointed out in ref. [S]. Concerning this issue, a new ingredient in this paper is 
SUSY breaking in systems in less than one-dimension. 
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into the two wells, and the value of the free energy is equal to that of the Gaussian SUSY 
matrix model. We think that it is deeply connected to the supersymmetric properties of 
the model, and further analysis will be reported in the next publication. 

On the other hand, we also consider the large- N and slowly varying limit for the same 
model. Then, solutions which exhibit spontaneous SUSY breaking are obtained. 

In sectional we consider the SUSY matrix model with a cubic W'(4>), where SUSY is 
expected to be preserved for finite N. It is found that spontaneous SUSY breaking takes 
place in the large- N limit with e — > even in this model. 

In section |6j we discuss on the SUSY matrix model with slowly varying Gaussian 
potential. Although the simple twist does not lead to SUSY breaking, deformations to 
^-dependent twists cause SUSY breaking in the large- N and slowly varying limit. 

Section [7] is devoted to summarize the results obtained so far and to discuss some 
future subjects. 

The path-integral expression of the partition function for SUSY quantum mechanics 
with twist a is derived in appendix El and appendix [B] is devoted to some fundamen- 
tal computation for the Gaussian SUSY matrix model. Finally, we give a proof of key 
equations concerning the SUSY restoration of the double-well matrix model at large- N 
in appendix O 

2 Preliminaries on SUSY quantum mechanics 

As a preparation to discuss large- N SUSY matrix models, in this section we present some 
preliminary results on SUSY quantum mechanics. 

Let us start with a system defined by the Euclidean (Wick-rotated) action: 



where is a real scalar field, ip, ip are fermions, and B is an auxiliary field. The dot means 
the derivative with respect to the Euclidean time t G [0,(3]. For a while, all the fields are 
supposed to obey the PBC. W(<f) is a real function of <p called superpotential, and the 
prime (') represents the 0-derivative. 

S is invariant under one-dimensional M = 2 SUSY transformations generated by Q 
and Q. They act on the fields as 



S 




) 



(2.1) 



Q<j) = Qip = 0, 
Qip = -iB, QB = 



(2.2) 



and 



Q<j> = -ip, 
-iB + 20, 



QB = 2ii, 



(2.3) 
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with satisfying the algebra 



Q 2 = Q 2 = 0, {Q,Q} = 2d t . (2.4) 
Note that S can be written as the Q- or QQ-exact form: 

S = Q J ' dt$^-B-(i + W\<j>))} (2.5) 

= QQ j dt fl^ + Wi^y (2.6) 

For demonstration, let us consider the case of the derivative of the superpotential 

W'(<P) = g^ + ^). (2.7) 

For fi 2 > 0, the classical minimum is given by the static configuration = 0, with its 
energy E = \g 2 ^ > implying spontaneous SUSY breaking. Then, B = —igfi 2 ^ from 
the equation of motion, leading to Qip, Qip ^ 0, which also means the SUSY breaking. 

For fi 2 < 0, the classical minima <p = ±a/— yU 2 are zero-energy configurations. It is 
known that the quantum tunneling (instantons) between the minima resolves the degen- 
eracy giving positive energy to the ground state. SUSY is broken also in this case. 

Next, let us consider quantum aspects of the SUSY breaking in this model. For later 
discussions on matrix models, it is desirable to observe SUSY breaking via the path 
integral formalism, that is, by seeing the expectation value of some field. We take (B) (or 
(B n ) (n = 1,2, • • •)) as such an order parameter. Whichever fi 2 is positive or negative, 
the SUSY is broken, so the ground state energy Eq is positive. Then, for each of the 
energy levels E n (0 < E < E x < E 2 <•••), the SUSY algebra^ 

{Q,Q} = 2E n , Q 2 = Q 2 = (2.8) 

leads to the SUSY multiplet formed by bosonic and fermionic states 

\b n ) = -±=Q\f n ), \f n ) = -±=Q\b n ). (2.9) 

As a convention, we assume that \b n ) and \f n ) have the fermion number charges F = 
and 1, respectively. Since the Q-transformation for B in (12.21) is expressed as 

[Q,B) = (2.10) 

in the operator formalism, we can see that 

(6 n |S|6„) = (f n \B\f n ) (2.11) 



2 In the operator formalism, Q, ip are regarded as hermitian conjugate to Q, tp, respectively. 
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holds for each n. Then, it turns out that the unnormalized expectation value of B 
vanishes^: 



(B)' = [ d(fields)5e- 

JpBC 



H\ IICJI !M JJ < 

'PBC 

Tr [B(-l) F e-? H ] 



J2(( b n\B\b n )-(fn\B\f n ))e'^ 



n=0 



= 0. (2.12) 

This observation shows that, in order to judge SUSY breaking from the expectation value 
of B, we should choose either of the SUSY broken ground states (|6 ) or |/ )) and see 
the expectation value with respect to the chosen ground state. It makes sense, because 
there is no transition between the two ground states via an arbitrary bosonic operator 
Ob (including B) from the conservation of the (—l) F - charge: 

(b \O B \f ) = 0. (2.13) 

The situation is somewhat analogous to the case of spontaneous breaking of ordinary 
(bosonic) symmetry in the sense that (12 . 1 3[) is reminiscent of the superselection rule. 

However, differently from the ordinary case, when SUSY is broken, the supersymmetric 
partition function vanishes: 



^PBC 



/ d(fields) e~ s = Tr [(-l) F e - pH ] (2.14) 

J PBC 

oo 

= £ (<W - </J/»>) e-'** 

n=0 

= 0, (2.15) 

where the normalization (b n \b n ) = (f n \f n ) = 1 was used. So, the expectation values 
normalized by Z PBC could be ill-defined [ID] . 

2.1 Twisted boundary condition 

To detect spontaneous breaking of ordinary symmetry, some external field is introduced 
so that the ground state degeneracy is resolved to specify a single broken ground state. 
The external field is turned off after taking the thermodynamic limit, then we can judge 
whether spontaneous symmetry breaking takes place or not, seeing the value of the corre- 
sponding order parameter. (For example, to detect the spontaneous magnetization in the 
Ising model, the external field is a magnetic field, and the corresponding order parameter 
is the expectation value of the spin operator.) 



3 Furthermore, (B n )' = (n = 1, 2, • • • ) can be shown. 



5 



We will do a similar thing also for the case of spontaneous SUSY breaking. For this 
purpose, let us change the boundary condition for the fermions to the TBC: 

iJ)(t + P) = e ia iJ)(t), ip{t + P) = e- ia iP(t), (2.16) 

then the twist a can be regarded as an external field. Other fields remain intact. As seen 
in appendix A, the partition function with the TBC corresponds to the expression (12.141) 
with (— 1) F replaced by (— e~ ta ) F : 

Z a = -e~ ia f d(fields) e~ s = Tr [{-e~ i,x ) F e~ pH ] (2.17) 

oo 

= J2{(bn\b n )-e~* a (f n \fn))e-^ 

n=0 

oo 

= (l-e~ ia )J2 e ^ En - ( 2 - 18 ) 

n=0 

Then, the normalized expectation value of B under the TBC becomes 
(B) a = ^Tr[B(-e-n F e-" H ] 

oo 

= Y S K 6 "^^) - e- ta (fn\B\f n )) e-P** 

a n=0 

agjj^ _ E^Q</n|B|/n>e-^ 



(2.19) 



Note that the factors (1 — e~ ia ) in the numerator and the denominator cancel each other, 
and thus (B) a does not depend on a even for finite ft. As a result, (B) a is equivalent to 
the expectation value taken over one of the ground states and its excitations {|& n )} (or 
{|/n)})- The normalized expectation value of B under the PBC was of the indefinite form 
0/0, which is now regularized by introducing the parameter a. The expression (12.191) is 
well-defined. 

On the other hand, from the Q-transformation xjj = [Q, 0], we have 

(U#n> = (fMfn) + -^=(fnWn)- (2.20) 

The second term is a transition between bosonic and fermionic states via the fermionic 
operator if), which does not vanish in general. Thus, differently from (B) a , the expectation 
value of <p becomes 

1 oo 
n=0 

EZ (fn\4>\fn)e-^ 1 EZoimK)^-^ 



■PE„ 



(2.21) 
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When {f n \if)\b n ) 7^ for some n, the second term is a-dependent and diverges as a — > 0. 
The divergence comes from the transition between \b n ) and |/ n ). Since the two states are 
transformed to each other by the (broken) SUSY transformation, we can say that they 
should belong to the separate superselection sectors, in analogy to spontaneous breaking 
of ordinary (bosonic) symmetry. Thus, the divergence of {4>) a as a — > implies that the 
superselection rule does not hold in the system Q 

2.2 Discretized SUSY quantum mechanics 

In this subsection, we consider a discretized system of (12. ip . namely the Euclidean time 
is discretized as t = 1, ■ ■ • , T. The action is written as 

S = Q^'W { ^B(t) - (^ + 1) - #) + (2.22) 

T 

= E L „ 

+ *(()W(t + l)-#)+r(#)Wt)} , (2.23) 



t=l 



~B(tf + iB(t) {cf>(t + 1) - cf>(t) + W'(4>(t))} 



where the Q-supersymmetry is of the same form as in (12. 2p . As is seen by the Q-exact 
form (I2.22p . the action is Q-invariant and the Q-supersymmetry is preserved upon the 
discretization [TT]. On the other hand, the Q-supersymmetry can no t be preserved by the 
discretization in the case of T > 2. 

When T is finite, the partition function or various correlators are expressed as a finite 
number of integrals with respect to field variables. So, at first sight, one might expect 
that spontaneous breaking of the SUSY could not take place, because of a small number 
of the degrees of freedom. In what follows, we will demonstrate that the expectation is 
not correct, and that the SUSY can be broken even in such a finite system. 

2.2.1 T = 1 case 

First, let us consider the simplest casdfl T = 1. As before, the boundary condition of the 
fermion is twisted with the phase a: 

0(2) = 0(1), ^(2) = e*V(l)- (2-25) 



4 For higher dimensional systems defined on finite spatial volume, such phenomena would also happen 
in general. 

5 For T = 1 with the PBC, the action is same as the dimensional reduction of (12. ip . Then, the 
Q-supersymmetry 

Q4> = Qi> = 0, 

Qtp = -iB, QB = (2.24) 

becomes the symmetry of the action again. 
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The partition function 



Z a = — [ dBd(pdijdije- Sa , (2.26) 

27T J 

S a = X -B 2 + iBW\(t))+^{e ia -1 + W"(4>))tp (2.27) 



with the superpotential (12.71) is computed to be 

1 



( e - _ i + W "^)) e -^'W 2 



= (e ia -l)C, (2.28) 
C = -— / d<Pe^ w '^ 2 . (2.29) 



C is positive definite, and it is same as in the continuum case that Z Q approaches to zero 
as a — ► 0. 

An analog to the unnormalized expectation value (12.121) with the PBC is 



(B)' = — [ dBd<pdi)di)Be~ Sa = 

d(j)W'((f))W"((f)) e -5 w 'W 2 



'2n J-oo 

However, the normalized expectation value (B) a becomes 



2tt 

• poo p\ 

/ ..a -iwrW = o. (2.30) 



(B) a = -t^TT- I dBd<f>d^d^Be- s - 

Z a ZlT J 

d<pW\<P) (e ia - 1 + W" ((/>)) e-^' w2 



1 -i 



oo 

POO 



/2tt y_ 

-i 1 



%rCJ_ 



oo 



#W(0)e-^'W 2 . (2.31) 



Again, the factors (e JQ — 1) from the numerator and the denominator cancel each other, 
and the result does not depend on a. For fi 2 > 0, (12.311) is not zero, implying SUSY 
breaking. At the classical level, B = —igfi 2 from the equation of motion, which coincides 
to ( 12.311) evaluated at the potential minimum = 0. This shows that the introduction 
of the external field a gives a correct regularization. Also, for fi 2 < 0, generically (12.311) 
does not vanish^, and it is found that the SUSY is spontaneously broken. 



6 Because 

B n = Q(i^B n - 1 ) (n = l,2,.-.), (2.32) 
we can conclude spontaneous SUSY breaking by showing (B n ) a ^ for some n. Even when (B) a = 
holds accidentally for /j 2 < 0, we could show (B n ) a ^ for some n and thus SUSY breaking. 
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The unnormalized expectation value with the TBC 

■ poo 

{B)' a = (e ia -i)—L d<\>W\<\>)e--* w '^ (2.33) 

vanishes in the limit a — > 0, and no singular behavior can be seen there. For the nor- 
malized expectation value, however, both of the numerator (B) and the denominator 
Z a approach to zero as a — ► 0. The dependence of a is canceled between the numera- 
tor and the denominator, and then the nonvanishing result arises. In other words, the 
external field a regularizes the indefinite form of the expectation value under the PBC: 
(B) =0/0, and leads to the nontrivial result. Thus, we conclude that, even in discrete 
systems defined by a finite number of integrals, SUSY can be spontaneously broken when 
the partition function vanishes. Note that it is totally different from the spontaneous 
breaking of ordinary (bosonic) symmetry, where the partition function never vanishes. In 
general, some non-analytic behavior is necessary for spontaneous symmetry breaking. For 
the SUSY breaking in the finite system, the non-analyticity is supplied by the vanishing 
partition function. 

As for the expectation value (<j>) a , similarly to the case (12.211) . we obtain 

{<t>) a = 7^ / dBd^di^dtPcPe-^ 
Z a lix J 

1 1 1 

e ia - 1 ^/2nC 

Since <pW"(<p) = 2g4> 2 , the integral of the right-hand side does not vanish, and (<p) a 
diverges as a — > 0. 

In what follows from section [31 we will consider various large- iV matrix models anal- 
ogous to the models presented so far. Interestingly, the singular behavior of (<f)) a as seen 
above does not appear there. It can be understood that the tunneling between separate 
broken vacua is suppressed by taking the large- N limit, and thus the superselection rule 
works. Note that the large-iV limit in the matrix models is analogous to the infinite 
volume limit or the thermodynamic limit of statistical systems. In fact, this will play an 
essential role for restoration of SUSY in the large-iV limit of the matrix model with a 
double-well potential. 




d<f)<J)W"(<J))e-^ w ' w2 . (2.34) 



General W'(<p) case Similar analyses can be done for W'(<j)) given as a general poly- 
nomial of the degree k: 

W'((j)) = g k <j) k + g k -i<p k - 1 + ■ ■ ■ + g . (2.35) 

Because 

^^(^■W^ ™ n (2,6, 

(from the Nicolai mapping |12j). and (12.331) holds for general W'(<f>), we see that the result 
of (B) a does not change from (12.311) for even k, while lim Q ^ {B) a = for odd k. (Note 
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that \im a ^ Z a ^ for odd k, since (I2.36P is nonzero.) Thus, for the discretized system 
with the superpotential W'(<f>) of the degree even (odd), the SUSY is broken (preserved), 
which is same as the conclusion for the continuous SUSY quantum mechanics [Hj. 

2.2.2 General T case 

It is straightforward to extend the above discussion to the case of general T. 
Expressing as S a the action ( 12. 23ft under the TBC 



the partition function 

Z a 

is computed to be 

Z a 
Cj 1 



0(T + 1) = 0(1), ^(T + l) = e i >(l), 



-i\ T r T 

— ) J Yl {dBif) d <t>(t) di/){t) #(*)) e ~ Sa 



t=i 



(-I f (1 - e ia ) a 

T 

n 



,t=i 



e 2- 



(2.37) 

(2.38) 

(2.39) 
(2.40) 



Here we used 



t=i 



r \TZ=xm+V-4>(t)+w'm))f = o (2.41) 



for the superpotential ( 12. 7ft . which is derived from the Nicolai mapping [12]. (Note the 



factor 



nli (-i + w"m))) 



is equal to the fermion determinant under the 



PBC.) Also, Ct is positive definite. 

Similarly, for the normalized expectation value 

rp rp 

(B(t))a = 4- ( ^ ) / fl #(*)) e" 



[2A2) 



we use the Nicolai mapping to have 

(B(t))a = 



H) (0(t + l)-0(t) + H/'(0(t))) 
xe -| £f =1 (0(m)-</>(t)+w" W*))) 2 



1 

Ct 



n 



d0(t) 



H) (0(t + i)-0(t) + H/'(0(t))) 
xe -|E?=i(^+i)-^)+W(0(t))) 2 _ 



(2.43) 
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The factor (— 1) T (1 — e ia ) was canceled, and (B(t)) a does not depend on a, again. The 
result (12.431) is finite and well-defined. By using the Nicolai mapping, it is straightforward 
to generalize this result to the case of (12.351) . We find that (12.431) holds and it is finite 
and well-defined for even k, and that lim Q ^o (B(t)) a = for odd k. 



No analog of Mermin-Wagner-Coleman theorem for SUSY As claimed in the 
Mermin- Wagner-Coleman theorem P, [7] , continuous bosonic symmetry cannot be spon- 
taneously broken at the quantum level in the dimensions of two or lower. In dimensions 
D < 2, although the symmetry might be broken at the classical level, in computing quan- 
tum corrections to a classical (nonzero) value of a corresponding order parameter, one 
encounters IR divergences from loops of a massless boson. It indicates that the conclu- 
sion of the symmetry breaking from the classical value is not reliable at the quantum level 
any more. It is a manifestation of the Mermin-Wagner-Coleman theorem. 

Here, we consider whether an analog of the Mermin-Wagner-Coleman theorem for 
SUSY holds or not. Naively, since loops of a massless fermion would be dangerous in 
the dimension one or lower, we might be tempted to expect that SUSY could not be 
spontaneously broken at the quantum level in the dimension of one or lower. However, 
this expectation is not correct. Because the twist a in our setting can also be regarded 
as an IR cutoff for the massless fermion, the finiteness of (12.431) shows that (B(t)) a is free 
from IR divergences and well-defined at the quantum level for less than one-dimension. 
(For one-dimensional case, (12.191) has no a-dependence, thus no IR divergences.) 

We can see it more explicitly in perturbative calculations. Let us consider the super- 
potential (12.71) with /i 2 > 0, where the classical configuration 0(t) = gives B{t) = —igfi 2 . 
If the theorem holds, quantum corrections should modify this classical value to zero, and 
there we should come across IR divergences owing to a massless fermion. Although we 
have obtained the finite result (I2.43p . the following perturbative analysis would clarify a 
role played by the massless fermion. We evaluate quantum corrections to the classical 
value of B{t) perturbatively. Under the mode expansions 

(T-l)/2 



0(t) = <Pn^ nt,T With 

V ^ n =-(T-l)/2 
(T-l)/2 

m = -4 £ ^ 2 ™+^, 



n=-(T-l)/2 
(T-l)/2 



m = -L Yl i> n e- t{2 ™ +a)t/T , (2.44) 



n=-(T-l)/2 



7 For theories in finite volume, the massless fermion might be regarded as "would-be Nambu-Goldstonc 
fermion" , because the Nambu-Goldstone fermion is defined for theories in the infinite volume and its 
concept is ill-defined for finite volume [5]. 
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free propagators are 



-nfm 



free 4 S in 2 (^)+M 2 ' 



with M 2 = 2g 2 fi 2 . Here we consider the case of odd T for simplicity of the mode expansion. 
Note that the boson is massive while the fermion is nearly massless regulated by a. Also, 
there are three kinds of interactions in S a (after B is integrated out): 



t=l 

T T 

v 3B = 9m 2 w + 1) - m) , v 3F = j2 2gmmm- (2.46) 
t=i t=i 

We perturbatively compute the second term of 

(B(t)) a = -ig» 2 - i (g^t) 2 + <t>{t + 1) - <t>{t)) a (2.47) 

up to the two-loop order, and directly see that the nearly massless fermion ("would-be 
Nambu-Goldstone fermion" ) does not contribute and gives no IR singularity. It is easy to 
see that the tadpole (<p(t + 1) — <p{t)) a vanishes from the momentum conservation. For 
—i (g<p(t) 2 ) a , the one-loop contribution comes from the diagram (IB) in Fig. [T], which 
consists only of a boson line independent of a. Also, the two-loop diagrams (2BBa), 
(2BBb), (2BBc) and (2BBd) do not contain fermion lines. The relevant diagrams for the 
IR divergence at the two- loop order are the last four (2FFa), (2FFb), (2BFa) and (2BFb), 
which are evaluated as 

(T-l)/2 / \ 2 



(2FFa) 

(2FFb) 
(2BFa) 




1 1 



e i(2irk+a)/T _ ]_ gi(27r(m+fc)+a)/T _ ]_ : 
2 



T 2 M 2 Z-J I 4sin 2 t™\ +M 2 / 4sin 2 /m) + M 2 

m=— IT— 1)/2 \ \ i / / \ i / 

1 



m=-(T-l)/2 
(T-l)/2 

X e i(2Trk+a)/T _ I 

fc=-(T-l)/2 



(2BFb) = £ 1- ' ^ 



T 2 M 4 ^ I 4 sin 2 (^) + M 2 ^ e i(2nk+ a )/T _ x 

m=-(T-l)/2 \ \ T / / fe=_(r-l)/2 




(2.48) 
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(IB) (2BBa) (2BBb) 




(2BFa) (2BFb) 



Figure 1: One- and two-loop diagrams. The crosses represent the insertion of the operator 
—ig4>{t) 2 . The solid lines with (without) arrows mean the fermion (boson) propagators. 
(IB) is the one- loop diagram, and the other eight are the two- loop diagrams. The diagrams 
with the name "FF" ("BB") are constructed by using the interaction vertices V^p twice 
(V4 once or V 3 b twice), and those with "BF" are by using each of V 3 b and V 3 p once. 
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Each diagram is singular as a — > due to the fermion zero-mode, however it is remarkable 
that the sum of them vanishes: 



(2FFa) + (2FFb) + (2BFa) + (2BFb) 



Ag 3 



T 2 M A 



1 



m=l 



M 2 



with 



F(m) 



E 

k=l 
T 

E 

fc=l 



1 + 



4 sin 2 (^) +M 2 



1 



3 j(27r(m+fc)+a)/T _ ]_ / e i(2irk+a) /T _ I 



e i(2-Kk+a)/T 



-i(2-Kk+a)/T 
_ e i2nm/T 



-i{2-nk+a)/T 
_ p -i2irm/T 



T 

e -i(27rfe+a)/T 

fe=l 



o. 



(2.49) 



(2.50) 



Thus, the two-loop contribution turns out to have no a-dependence, and the quantum 
corrections come only from the boson loops which are IR finite, that is consistent with 
(12.43j) . Since the classical value —igfJ? = is regarded as 0(g~ 1 ), and £-loop con- 

tributions are of the order 0{g 21 ^ 1 ), the quantum corrections can not be comparable to 
the classical value in the perturbation theory. Thus, the conclusion of the SUSY breaking 
based on the classical value continues to be correct even at the quantum level. 



3 Overview of SUSY matrix models with TBC 

In this section, we define a matrix model analog of the discretized SUSY quantum me- 
chanics with the TBC as discussed in the previous section, and give an overview of its 
general features in the large- iV limit. 

Let us begin with a matrix-model analog of (12.221) 

T 

S = Qj2 mT 

t=i 

with the Q-supersymmetry nilpotent (Q 2 = 0): 

Q<p(t) = m, q#) = o, 

Q$(t) = -iB(t), QB(t) = 0. (3.2) 

All the variables are N x N hermitian matrices and the derivative of the superpotential 
is a polynomial of </> given as 

wm)) = Y,g P m p - (3.3) 

v 



1) - <j>{t) + W\<t>{t))) 



(3.1) 
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Then, the action is 



s = j2 Nti 



t=i 



l -B(t) 2 + iB{t) {<f>{t + 1) - cf>{t) + W'{<t>{t))} 



p k=0 



(3.4) 

We will focus on the simplest case T = 1 with the TBC (I2.25P again 0. Under the action 



S n = Ntr 



-B 2 + iBW'{<t>) + i> (e ia -1)^ + J29pJ2 ^<f~ l 

p k=0 



(3.6) 



the partition function is defined by 



Z a = (-1)^ J d N2 Bd N2 <f)d N2 ijd N2 ije- s «. 
For convenience, we fix the normalization of the measure as 



(3.7) 



d N^ e -Ntr(^) 



d N ^d N ^e 



-TVtr (ipip) 



(3-8) 

In order to construct SUSY matrix models which exhibit spontaneous SUSY breaking 
in the large- N limit, we introduce a small positive constant e and consider the situation 
that the derivative of the superpotential W'(</>) depends on cf> only through the combination 
of e</>, i.e. 

(3.9) 



W\ct>) = h{e < P)=Y J 9p^ 



with g p constants of order 1. Later we will discuss spontaneous SUSY breaking in the 
case that e tends to in the large- N limit. (For example, e scales as 0(N~ U ) with u > 0.) 
Since e~ l corresponds to the length scale of in which W'(<f>) varies by an amount of 
0(1), the potential W'((f)) becomes more and more flatter or slowly varying as N — > oo@. 



8 Similarly to section [2721 for T = 1 with the PBC, the action is supersymmetric under Q as well as 
under Q for a general superpotential, and it can be expressed as 



S' = ()()w [ ^ + W(cf>) 



(3.5) 



The Q, Q transformations are of the same form as in (|2.2j) and (I2.24|) . and nilpotent (Q 2 = Q 2 = {Q, Q} = 
0). 

9 In the case of large T, we can consider the continuum limit where the lattice spacing of the discretized 
time a — > with aT fixed. It would be interesting if e can be also related to a in the large-iV limit. It 
is similar to the situation of the double scaling limit in the matrix models for two-dimensional quantum 
gravity [TJ [HI US] . 



15 



Instead of taking the limits iV — > oo and e — > simultaneously, we can consider e as 
independent of N, and can take the limits separately. Note that the external field a is 
always turned off at the final step because it is for detecting spontaneous SUSY breaking 
of the system obtained after the large- N limit and the e — > limit. In our model, it is 
easy to check that irrespective of the order of the limits N — > oo and e —>■ 0, we will have 
the same large- N saddle point equation at the final step a — > 0, and thus we obtain the 
same result. Therefore, we do not have to care about the order of these limits. 

After integrating ip, ijj, B and rescaling as <p ~~ > e_ V) t ne partition function becomes 



Z n 



z' = 



eN 2 n ai 



]N 2 



det 



(^-lJlOl + eJ^Z)' 



fc=0 



-JVtrfih^) 2 ! 



(3.10) 



The expectation values of -^tr B r 



[n 



1,2, 



are expressed as 



Z n 



,JV 2 



x det 



( e »*-l)l®l + eJ> P E 



fc=0 



-JVtr[ifc(^) 3 l 



(3.11) 



To figure out how the large-iV and slowly varying limit causes spontaneous SUSY 
breaking, we first regard e as independent of N, and consider the case of large N but finite 
e in section 13.11 Then, the large- N and slowly varying limit is discussed in sections 13.21 
and I 



3.1 Large- N solutions for finite e 



We consider large- N solutions for the system (13.101) . (13.111) with e kept finite. Here, 
we do not take the e — > limit, and the external field a is turned off after the large- 
N limit. (13.101) can be expressed as integrals with respect to the eigenvalues of (A» 
(i = l,...,JV)): 



Z' a = C N [ (f[d\A A(A) 2 fl 

J \i=l J i,7=l 



p-1 



k=0 



(3.12) 
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with a numerical constant and A(A) = n*>j — ^j)- m the large-^ limit, the 
integrals can be evaluated at the saddle point configuration as 

Z' a = C N exp[-N 2 (F a + 0(N- 2 ))], 
F a = — J dx dyp{x)p{y) P In \x — y\ 

- I dx dyp{x)p{y) In (e ia - 1 + e £ g p ^ x k y p ~ l - k 

J \ p k=0 / 



j dx-h{x) 2 p{x) 



(3.13) 



with Cn = exp [|iV 2 + (9(iV )] obtained in (IB.8P in appendix [Bl The eigenvalue distri- 
bution p{x) = jj Eili <K X ~ ^i) satisfies the large- N saddle point equation: 



dyp{y)P- 



P- 1 u^k-U.v-l-k 



+ / dyp{y) 



'■P- 1 rr-k..V-1-k 



x ~y J ^-i + ej2 p g P T: p k=o^y p - 

Finally, after the limit a — > 0, we obtain 



-h(x)ti(x). (3.14) 



dyp{y)Y- 



x-y 



+ / dyp(y) 



p- 1 „jfe„p-i-* 



-h(x)ti(x). (3.15) 



Note that the e-dependence disappears in the equation. 

In terms of p(x) which solves (13. 15ft . the expectation values of -^tr B n {n = 1,2, 
are given by 

lim | lim / — tr B 



i / dxp(x)h(x) 



lim I lim ( — tr B 2 

a^O ViV->oo \ N 



lim I lim ( — tri? 3 

a^O V N-+oc \ N 



lim I lim ( — tr B 4 



dxp(x)h(x) 



3z / dx p{x)h(x) + i / dx p(x)h(x 



2 - 4y dxp(x)h(x) 2 - 2 
+ / dx p(x)h(x) 4 , 



-i 2 



da; p(x)h(x) 



\ . (3.16) 

In section HH1 we explicitly solve the saddle point equation for the quadratic W((f>): 
W'((f>) = h(e<f)) = (e<f)) 2 + p? by the technique of the 0(n) model on a random surface 

with ti = -2 [na oil rig. 
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Here it is important to recall that, when both of e and N are finite, basically we will 
obtain the same result as in the discretized SUSY quantum mechanics. Namely, when 
W(<f>) is a polynomial of odd (even) degree, the SUSY is preserved (broken). However, 
a new feature can arise in the large- N limit even if e is kept finite. It is analogous to 
the infinite volume limit of SUSY field theory [5]. There, in the case of the double- well 
potential, the SUSY breaking is triggered by instanton effects in finite volume, however 
the instantons are suppressed in the large volume limit, and the SUSY becomes restored. 
We will find in section |4~T1 that the large- N limit plays the same role as the infinite volume 
limit and restores the SUSY for the same potential. 



3.2 Large- N SUSY breaking under slowly varying potential 



In the large-iV and slowly varying limit, N — > oo together with e — ► 0, we can discard 
the e-dependent terms in the determinant factor in (13.101) . Then, the system is essentially 
reduced to the ordinary one-matrix model with respect to <fi. In terms of the eigenvalue 
distribution p(x), the partition function Z' a is given in the large- N limit as 



Z' a 
Fa 



(e ia - if 2 C N exp [-N 2 (F Q + OiN' 2 ))} 



dx dyp(x)p(y) Pin \x — y\ + / dx-h(x) 2 p(x) 



with p(x) satisfying the saddle point equation 

1 



dyp(y)P- 



x-y 



-h(x)h'(x) 



(3.17) 



(3.18) 



Note that the equation does not depend on a. It can be solved in the standard method 
[T§] by assuming a finite support for p(x). For example, given such a solution p(x) with 
the support x G [a, b], (13.171) can be expressed by integrating (13. 18j) as 



dxp{x) In \b — x\ + I dx-h(x) 2 p(x) + -h(b) 2 . 



(3.19) 



As for the expectation values of j^ti B n (n = 1, 2, 
1 



we obtain 



-MB 



N 



F" 4 



dxp(x)h(x), 
dxp(x)h(x) 2 , 



3i / dx p(x)h(x) + i / dx p(x)h(x 



2 — 4y dxp(x)h(x) 2 - 2 
+ / dx p(x)h(x) 4 , 



dx p(x)h(x) 



(3.20) 
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which do not depend on a similarly to the situation seen in the (discretized) SUSY 
quantum mechanics in section [21 
Also, 

^tr(( e 0)")^ =Jdxx n p(x) (n=l,2,...) (3.21) 

do not depend on a and exhibit no singular behavior as a — > 0. It is different from 
the situation of the (discretized) quantum mechanics, where (<f>) a brows up in the limit 
a — > as in (I2.2ip and (I2.34p . It implies that the large- N limit suppresses the quantum 
tunneling to make the superselection rule work in the matrix model. Although (^tr (0 n )) Q 
grows as e~ n in the slowly varying limit, any singular behavior does not appear for the 
"renormalized" variable e0. 

Now it is clear from (13.101) and (13.111) how SUSY breaking happens in the e —>■ limit. 
Since the large- N and slowly varying limit amounts to dropping the fermion interactions, 
there is generically no reason to expect that the SUSY persists after turning off a. It 
will be so regardless of the form of W'(<p), except the Gaussian case W((fi) = h(ecf)) = ecj) 
where the fermion determinant is constant. In the following sections, for some concrete 
examples of W'((f>), we compute the expectation values of -^tr B n {n = 1,2, •••) from 
which we observe spontaneous SUSY breaking, and see that this is indeed the case. In 
particular, it is quite interesting that SUSY is broken even in the case of W'(<f)) with odd 
degree, for which the SUSY is preserved in the (discretized) SUSY quantum mechanics. 
Since the finite N case will lead to the same conclusion as in the discretized quantum 
mechanics, our model provides some new examples of spontaneous SUSY breaking in 
that the a — > limit and the large- N limit do not commute. It is more analogous to 
spontaneous breaking of ordinary bosonic symmetry. 

As seen from (13.101) and (13. lip , without the twist a, we cannot take the e — > limit in a 
well-defined manner. In a sense, this is as usual: in order to discuss spontaneous symmetry 
breaking, introduction of the external field is essential, because it lifts degeneracy of vacua 
and makes the notion of the vacuum expectation value well-defined as stressed repeatedly 
in the previous section. 



3.3 Another deformation 

So far, we introduced the external field a to see the spontaneous SUSY breaking. It can 
be regarded as a deformation of the model with inserting the operator, which changes the 
boundary condition from the PBC to the TBC by the twist a. The deformation slightly 
breaks the SUSY of the original model. Then, what we saw in section 13.21 is interpreted 
as a statement that, if the limit a — > is taken after the large- N and slowly varying limit, 
the SUSY breaking effect remains even after turning off the deformation a. 

To see how the model is sensitive to the deformation, we here consider another de- 
formation corresponding to a specific twisted boundary condition, where the twist also 
depends on the variable 0: 

0(2) = 0(1), i/j(2) = e^trv^i))^^ (3 >22 ) 



19 



where V(<j>) is a function of <fi. Then, the action reads 



S n = Ntl 



p-1 



-B 2 + iBW'((f>) + iJ (j a > v (e<t>) + ^ e p g p 



k=0 



. (3.23) 



Note that the twist depends on <fi through the combination e<p, which is the same as 
the ^-dependence of the potential W'(<f>) = h(e(f>). We consider such a setting in order 
to make well-defined the limit a — > after the large- iV and slowly varying limit. The 
^-dependent twist could also be interpreted as some operator inserted at the boundary. 

Similarly, after the rescaling — ► e _1 0, the partition function and the expectation 
values of -^tr B n (n = 1, 2, • • • ) become 



Z' = I d N2 6 det 



giaitry(^) 



p-1 
p k=0 



up-l-k 



xe 



-Ntr \\h{^) 2 



(3.24) 



N 



-tr B' 



1 



d N B d 



x det 



p-i 

(e ia ^ y W - l) 1 ® 1 + e 9 P 



ip-i-k 



xe 



-7Vtr[±/i(«« 2 l 



(3.25) 



In the large- N and slowly varying limit, the e-dependent part in the determinant factor 
can be neglected. Then, the partition function Z' a at a large- N saddle point is expressed 
as 



Z' a = C N exp[-N 2 (F a + 0(N- 2 ))}, 

F a = - J dxdyp{x)p{y)Y\n\x-y\ -\ n (e ia $ dxV{x)p{x) - l) + Jdxhi{x) 2 p{ 
where the eigenvalue distribution p(x) is determined by 



dyp(y)P 



1 hiaV'(x) 1 



+ 



x — y 1 — e 



h(x)h'(x) = 

iac 2 



X), 

(3.26) 
(3.27) 



with c = J dyV(y)p(y). Note that the second term remains nonzero and contributes as a 
part of the potential lj even after the limit a — > 0: 



dyp(2/)P— ^— = ^(x)/i'(x) - — V'(x). 

x — y 2 2c 



(3.28) 



in 



This would be a peculiarity of our model with T = 1. 
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Therefore, our problem is reduced to that of the one-matrix model. In this case, (13.281) 
should be solved consistently with c = J dyV \y) p(y) . In the following sections, for various 
W'((j)) and V(e(j)), we present the solutions. As a relevant quantity for the planar free 
energy after the limit a — * 0, we will consider 

F = lim (F a + ln(ia)) 

a— >0 

= — In c ~ J dxdyp(x)p(y) Pin \x — y\ + J dx -h(x) 2 p(x) . (3.29) 

Given such a solution p(x), for example with a support x G [a, b], we can express (13.291) 
as 

F = -hie- J dxp(x)ln\b- x\ + ^(c-V{b))+ dx p{x)\h{xf + \h{bf . (3.30) 

The expectation values (^tr B n } {n = 1, 2, • • • ) after a turned off are expressed in the 
same form as (I3.20p . Note that it is finite and well-defined in the limit a — > 0. For the 
discretized quantum mechanics with the slowly varying potential W '((/)) = h(e<f)) of odd 
degree and the 0-dependent twists aV(e4>), it will be observed that some of (B n ) a are 
singular as a — > after the slowly varying limit e — > 0. Since this kind of pathological 
behavior does not appear in the case of the large-iV matrix model, it shows further 
importance of the large- iV limit to construct the models. 

For the same reason as in section 13.2} we can say that spontaneous SUSY breaking 
happens for any W'((p) in the case of the ^-dependent twist (I3.22j) . Since, differently 
from the constant twist (12. 25ft . the twist V(<p) now takes part in the saddle point equa- 
tion (I3.28p . the twist V(e<j)) causes spontaneous SUSY breaking even for the Gaussian 
case W'((j)) = h(e4>) = ecfi, as we will see in section El 



4 Large- N </> 4 SUSY matrix model 

In this section, we concretely examine the case with 

W\<P)=h(e<P) = (e<P) 2 + p 2 (4.1) 

in the setting of section [3l 



4.1 Large- N solutions for finite e 

Here, we will see that the system (I3.10p with (14. ip for general finite e has large- N saddle 
points, which is expected to preserve SUSY for (some region of) p? < 0. For finite N, 
on the other hand, the system will take place spontaneous SUSY breaking, as we saw 
explicitly for the N — 1 case in section 12.21 Thus, this presents an example of matrix 
models, where the SUSY is broken for finite iV but restores in the large- iV limifl 

11 Such an example for a large- N vector model has been presented in Ref. [1]. 
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The partition function reads 



7 - — 7' 



f / N \ N 

Z' a = C N (HdxA A(A) 2 J] [e ia - 1 + e(Aj + Xj)] e^EMC** 3 ) , (4 . 2 ) 
J \i=i J i,j=i 



and the saddle point equation ( 13. 15ft becomes 

dyp(y)P^— + I dy = x 3 + p 2 x. (4.3) 

x - y J x + y 

Let us consider the case p 2 < 0, where the shape of the potential is a double- well 
4.1.1 One-cut solution 

First, we find a solution corresponding to all the eigenvalues located around one of the 
minima A = +a/— p 2 - Assuming the support of p(x) as x G [a, b] with < a < b, the 
equation (14.3ft is valid for x G [a, 6]. 

Following [15] . we introduce a holomorphic function 

F(z)= [ dy-^-, (4.4) 

A 2-2/ 

which satisfies the following properties: 

1. F(z) is analytic in z G C except the cut [a, b] . 

2. is real onzeR outside the cut. 

3. For 2 ~ oo, 

F(z) = i + 0(i). 

4. For x G [a, 6], 

F(x ± «0) = F(—x) + x 3 + p 2 x =F mp(x). 

Note that, if we consider the combination 



F-{z) = \{F{z)-F{-z)), (4.5) 

then the properties of F_(z) are 

1. F-(z) is analytic in z 6 C except the two cuts [a,b] and [—6, —a]. 

2. is odd (F_(—z) = —F-(z)), and real onzeR outside the cuts. 
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3. For z ~ oo, 

F_(*) = l + 0(£). 

4. For x G [a, b], 

F_{x ± zO) = ~ (x 3 + /i 2 x) =F i\p{x). 

These properties are sufficient to fix the form of F-(z) as 

F_(z) = ~ (z 3 + fi 2 z) - l - z ^( z 2 -a 2 ){z 2 -b 2 ) (4.6) 

with 

a 2 = -2 - /i 2 , 6 2 = 2 - /i 2 . (4.7) 

Since a 2 should be positive, the solution is valid for — ji 2 > 2. The eigenvalue distribution 
is obtained as 

p(x) = -\/{x 2 -a 2 )(b 2 -x 2 ). (4.8) 

Applying (14.81) to (13.161) . it is seen that all the expectation values of -^tr B n {n = 
1,2,3,4) vanish. Moreover, we can show 

lim ( lim ( — txB n ) ) = (n = l,2,---) (4.9) 



by an inductive argument as in appendix O 

Also, the partition function can be calculated, and we see that it coincides with that 
of the Gaussian SUSY matrix model with g% > obtained in appendix [Bj 

lim ( lim -^-lnZ a ) = lim -L ln ( z c\ n >n) = 0- (4.10) 
These evidences convince us that the SUSY is restored at infinite N. 



4.1.2 Two-cut solutions 

Let us consider configurations that v + N eigenvalues are located around one minimum 
A = +\/—fJ> 2 of the double-well, and the remaining v_N(= iV — v+N) eigenvalues are 
around the other minimum A = — a/— /i 2 . Since the fermion determinant in the partition 
function (14. 2 j) is 

N N 

J] [e ia - 1 + e(Ai + A,)] = J] [ e " " 1 + 2eA *] II [ e " ~ 1 + < Xi + ^ > ( 4U ) 

i,j=l i=l «>J 

the configurations will contribute (— l) 1 '-^ to the sign of the partition function for small 
a. Thus, the sign is stable against variations around the configurations with v + or V- 
fixed. Therefore, in the large-iV limit the tunneling is suppressed and a smooth saddle 
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point is expected to be found. Because such a saddle point will form a two-cut solution 
at large N, let us investigate two-cut solutions. 

First, we focus on the ^-symmetric two-cut solution with u + = v- = ~, where the 
eigenvalue distribution is supposed to have a ^-symmetric support Q = [—b, —a] U [a, b], 
and p(—x) = p(x). The equation (14. 3B is valid for x G Q. Due to the Z 2 symmetry, the 
holomorphic function F(z) = J n dyj^ has the same properties as F_(z) in section 14.1.11 
except the property 4, which is now changed to 

F(x ± zO) = - (x 3 + fi 2 x) =F inp(x) for igSI. (4.12) 



The solution is given by 



F(z) = l - (z 3 + p?z) - \zyj(z* - a 2 ){z 2 - b 2 ), (4.13) 

p(x) = ^-\x\y/(x 2 -a 2 )(b 2 -x 2 ), (4.14) 

where a, b coincide with the values of the one-cut solution ( 14.71) . It is easy to see that, 
concerning Z 2 -symmetric observables, the expectation values are same as the expectation 
values evaluated under the one-cut solution. In particular, the proof in appendix O can 
be applied also here, to obtain 

lim ( lim ( — tr B n ) ) = (n=l,2,---)- (4.15) 



a^O \N->oo \ 

The partition function with the sign factor dropped is the same as that evaluated by the 
one-cut solution: lim a ^ (limjv^oc, In |Z a |) = 0. 

It is somewhat surprising that the end points of the cut a, b and the absolute value 
of the partition function coincide with those for the one-cut solution, which is recognized 
as a new interesting feature of the supersymmetric models and can be never seen in the 
case of bosonic double-well matrix models. In bosonic double-well matrix models, the free 
energy of the Z 2 -symmetric two-cut solution is lower than that of the one-cut solution, 
and the endpoints of the cuts are different between the two solutions [201 12"T] . 

Next, let us consider general ^-asymmetric two-cut solutions (i.e., general v±). We 
can check that the following solution gives a large- iV saddle point: 

The eigenvalue distribution p(x) has the cut Q = [—6, —a] U [a, b] with a, b given by (14. 7ft : 

/ \ f — x J (x 2 — a 2 )(b 2 — x 2 ) (a<x<b) 
p{x) = { y y rr-, ?T773 k L . J ^ (4-16) 



\x 



a/(x 2 - a 2 ){b 2 - x 2 ) (—b<x<-a). 



This is a general supersymmetric solution including the one-cut and -^-symmetric two-cut 
solutions. The expectation values of Z 2 -even observables under this saddle point coincide 
those under the one-cut solution, and lim^o (linijv^oo 

again. Also, for the partition function, limc^o (liniA^oo ^2 In \Z a \) = 0. Thus, we can 
conclude that the SUSY matrix model with the double-well potential has an infinitely 
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many degenerate supersymmetric saddle points parametrized by at large N for 

the case p 2 < —2. It is totally different from the N — 1 case discussed in section 12.2. 1} 
where the partition function with the PBC vanishes and the SUSY is broken. 

For the general two-cut solution with (u + ,v-), the partition function appears to be 
expressed as 

lim ( lim — ln|Z a | J = lim -r=\n\Z Gv+ Z G j,_\ =0. (4.17) 

Here, Zq v± are the partition functions of the Gaussian SUSY matrix models with the 
matrix size v±N, whose superpotential is obtained by the Gaussian approximation of the 
original superpotential W'(x) = h(ex) = (ex) 2 + p 2 around its minima x = ±-a/— p 2 , 
respectively. Furthermore, since Z a and Zg, v+ Zg, v ^ have the same sign factor (— l)^ - ^, Z a 
can be evaluated by the product of the partition functions of the Gaussian SUSY matrix 
models Zq,v + Zg,v_ including the sign factor at large N. It is reasonable to attribute these 
remarkable features for general v± to the restoration of the SUSY. Details of this aspect 
will be reported in the next publication. 



4.2 Large-iV SUSY breaking under slowly varying potential 

Here, we show that the model with the same superpotential fl4.ll) exhibits spontaneous 
SUSY breaking in the large- N and slowly varying limit. 
The saddle point equation (13.181) becomes 

dyp(y)P—^— = x 3 + fi 2 x. (4.18) 
x-y 

One-cut solution Assuming the support of p(x) as x G [—a, a] and using the standard 
method [19j, it can be solved as 



p(x) -- 

7T 

2 2 

a — - 
3 



l - (x 2 + p 2 + v 7 ^^ 2 , (4.19) 
(-p? + ^/pJT6) . (4.20) 



From the requirement p(x) > 0, this solution is valid for p? > —\/2. 

Then, the expectation value of -^tr B becomes independent of a at large N: 



—trBj = —il dx (x 2 + p 2 )p(x) 

^ la J -a 



i\a 2 (a 2 + 2p 2 ) 2 . (4.21) 



8 

This is nonzero for p 2 > -y/2, meaning the SUSY is broken. It accidentally vanishes at 
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p 2 = — v2, but since 

1 



ca 

^tr5 2 ) =l~l dx (x 2 + p 2 ) 2 p(x) 



—a 



= 2^° for A* 2 = -V2, (4.22) 

the SUSY breaking is seen also at p 2 = —y/2. 

As a quantity relevant to the planar free energy (after subtracted by the divergent 
constant as a — > 0), we compute (13.171) : 

F = — J dx dy p{x)p{y) P In \x — y\ + J dx - (x 2 + p 2 ) 2 p{x). (4.23) 

Plugging the one-cut solution f)4.20p into (I3.19p . we obtain 

^o|„ t = -ln^ + | + ^ 4 +f^ 2 + ^. (4.24) 



2 8 36 r V24 72 

Two-cut solution For p 2 < —a/2, the solution develops to have two cuts. Among 
various two-cut solutions, the most relevant is a Z 2 -symmetric one [201 EI], where the 
eigenvalue distribution has two supports [—6, —a] and [a, b). The explicit form is 

p(x) = - \x\^(x 2 -a 2 )(6 2 -x 2 ) (4.25) 

7T 

with a 2 = -p 2 -V2,b 2 = -p 2 + y/2. 

Interestingly, the free energy for this solution is turned out to be /i 2 -independent: 

^o| two _ cut = ^ + ^ln2, (4.26) 

and the expectation values of B are 

itr B ) =0, (ItrflA =1. (4.27) 

The SUSY is broken also at the two-cut phase. The property of the free energy might be 
interpreted as some remnant of supersymmetric properties of the model. 

Phase transition Comparing -Pol one-cut an< ^ I two-cut a ^ ^ ne P nase boundary p 2 = 
—a/2, from ( 14.24p and (14.261) it is seen that F , j^Fq and d ^ 2 ^ Fo is continuous, but 

d 3 1 d 3 

lim -t^-^o = -=, lim -rr^ F Q = °> ( 4 - 28 ) 

M 2_»_V2+o d(/i 2 ) 3 At ^_ v ^_ rf(/i 2 ) 3 

implying that the phase transition is of the third order. 
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4.2.1 Another deformation 



We consider the same model but with the twist a changed to a 0-dependent function. As 
the function V(e<f>) in (13.221) . we choose the simplest ^-symmetric one V(e<f) = (ec/>) 2 for 
a technical reason. 

After the limit a — > 0, the corresponding saddle point equation (13.281) is 



dyp(y)P—^— = x 3 + ( p 2 - — ] x (4.29) 
x-y V c 2/ 



with c 2 = Jdyy 2 p(y). 



One-cut solution To solve the equation (14.291) . we assume the support of p(x) as 
x G [—a, a] from the Z 2 symmetry. The holomorphic function F(z) = J^ a dy is 



determined as 



1 \ / o 8 — a 



F{z) = z 3 +(p 2 --)z-(z 2 + ) Vz^a* (4.30) 



with 



Then, 



a 2 , 4 



° 2 = 32 (a +8) - (431) 



1 / n 8 — a 



7T 



p( x ) = - ( x 2 + ) v^ 2 ^ 2 , (4.32) 



and a is given as a solution of the equation 

3a 8 + 4/i 2 a 6 + 16a 4 + 32p 2 a 2 - 192 = 0. (4.33) 

Although the explicit form of a is somewhat complicated, it can be seen that < a 2 < 2V2 
when p 2 > —l/y/2. Then, p(x) > is satisfied, and the solution is valid. SUSY breaking 
is observed from the expectation value 

i^o (& (jj trB ) a ) = ~ l ^ + ^ * °' (434) 
In particular, when p = 0, we obtain a simple expression of the solution as 



a 4 = ?(VT6-l), c 2 = ^\M^-i(^o + 2 ). 

As a relevant quantity to the planar free energy, we compute (13. 29ft : 

F = — lnc2 — J dx dyp(x)p(y) Pin \x — y\ + J dx - (x 2 + w 2 ) 2 p(x). (4.36) 
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For the one-cut solution, we obtain 







-1/V5 - 


1 ln2+ 5 
4 ln2+ 8 


d F 

At , 2\- t ° 


y? = 


v 


1 

/7T' 

V2 


d " F 
d(/i 2 ) 2 


M 2 = 


-1/V2 


1 

3' 


d3 F 






8^ 


d(^) 3 


M 2 = 


-1/V2 


27 ' 



(4.37) 



Two-cut solution For /i 2 < — 1/ a/2, the relevant large- N solution has the Z 2 -symmetric 
two cuts Q = [—b, —a] U [a, b] (0 < a < b). The solution is 



c 2 = 2 1 -V 



vV + 4 
p(x) = — \x\\J {x 2 — a 2 )(b 2 — x 2 ). 



a = c 2 



V2, 



71 



b 2 = c 2 + y/2, (4.3£ 



(4.39) 



SUSY breaking is seen from 



lim lim ( — tr B 



(/i 2 + a/a^+4) o. 



(4.40) 



d 2 



Computing the free energy F for the two-cut solution, we find that F , ^jj-Pb an d 
-F at /j 2 = -I /y/2 coincide with (OTjl . but that 



d 3 



2\3' 



4a/2 



(4.41) 



is different. Thus, the phase transition at /i 2 = — l/v2 is of the third order. 



Remarks Comparing the result here with the one obtained for the constant twist, we 
see that the change of the deformation, namely the constant twist a changed to the 
^-dependent function a;4tr ((e0) 2 ), modifies the precise form of the one- and two-cut 
solutions for p(x), but their qualitative behavior and the order of the phase transition 
remain unchanged. It is considered that, because in the slowly varying limit the external 
field is dominant in the fermion bilinear terms, the model will have some sensitivity to 
the form of the twist even after it is turned off. 
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5 Large- N 6 SUSY matrix model under slowly vary- 
ing potential 

In this section, we consider the SUSY matrix model with the derivative of the superpo- 
tential 

W'(<j>) = h{e<j>) = (ex) 3 + p 2 (ex). (5.1) 

For finite N and finite e, it can be seen that SUSY is not broken, for example by using the 
Nicolai mapping. We show that even in this system, at the large- N and slowly varying 
limit (e — ► 0), the SUSY is spontaneously broken. 

For the constant external field a, the corresponding large- iV saddle point equation 
(13. 18ft becomes 

1 1 

dypiy) P = -(x 3 + p 2 x)(3x 2 + p 2 ). (5.2) 

x — y 2 

The one-cut solution with the support x G [—a, a] is obtained as 

p(x) = — (3x 4 + A 2 x 2 + Ac) Va 2 - x 2 , (5.3) 

where 

A 2 = V + la 2 , A =fi 4 + 2fi 2 a 2 + la 4 , (5.4) 
2 o 

and a solves the equation 

15a 6 + 24/i 2 a 4 + 8/iV - 32 = 0. (5.5) 
Then, ^tr5) vanishes due to the Z 2 symmetry, but 

jt^B 2 \ =1-/ dx (x 3 + u 2 x) 2 p(x) (5.6) 

la J —a 

is independent of a, and nonzero in general implying spontaneous SUSY breaking. This 
solution is valid for the case fx 2 > — 3 2 / 3 . At the boundary /i 2 = — 3 2//3 , the eigenvalue 
distribution takes the form 

p(x) = ^-(x 2 - 3- 1/3 ) 2 Va 2 - x 2 (5.7) 

with a = 2 ■ 3~ 1//6 . It touches the zero at x — ±3 -1 / 6 , and then 1/3 of the eigenvalues are 
distributed in each of the three regions [-2 • 3~ 1/6 , -3~ 1/6 ], [-3~ 1/6 , 3" 1/6 ] and [3~ 1/6 ,2 ■ 
3-1/6]. For p? < — 3 2 ^ 3 , the solution develops to have three cuts. 



p = case For the case p = 0, we obtain a simple expression: 

|) 1/3 , A 2 = 3 -a 2 , A = 9 -a\ (5.8) 
and 

1 \ f a 2 

— tiB 2 ) =1-/ dxx 6 p(x) = - ^ 0. (5.9) 

/ ex. J Q> ^ 
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5.1 Another deformation 



Similarly to section 14.2.11 we consider the deformation of the twist a corresponding to 
V(&p) = (e0) 2 in (I3T22D . 

From (13.281) the saddle point equation after a — > is 



dyp{y)P 



x — y 2 



li 



C-2 



x + 4/i 2 x 3 + 3x 



(5.10) 



with C2 = f dyy 2 p(y). Assuming the support of p(x) as x G [—a, a], we solve it as 

1 



= 7T" (3x 4 + A 2 x 2 + Aq) Va 2 - x 2 , 
Ztc 



where 



V + -a 2 , 



4 3 



A 2 

c 2 = 7^77 (15a 8 + 16p 2 a 6 + 64a 2 ) . 
25b 



— —a — pa 



And a is determined by 

(15a 6 + 24/i 2 a 4 + 8/i 4 a 2 - 32) (15a 6 + 16/iV + 64) = 16 • 256. 

p = case For p = 0, eq. (I5.13P is easily solved to obtain 



64 
15 



1/3 



a ° 2 

Ao = -a 2 , 



An 



4 



2~ 1 — ha 2 ' 
Then, ^trB n ^ a with n odd vanish from the Z 2 symmetry, but 
1 



c 2 = -a 



lim lim ( —trB 



\N^oo \ N 

lim | lim / — tr5 4 

V^^oo \ ]V 



dx p(x)x = 0, 



2 — 4 / dx p(x)a; 6 — 2 I / dxp(x)x 3 ] + / dx p(x)x 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



,12 



34 
225 



^0, 



(5.15) 



implying spontaneous SUSY breaking. 
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6 Large- N Gaussian SUSY matrix model under slowly 
varying potential 

In the previous sections, we have investigated the SUSY breaking in the large-iV and 
slowly varying limit for W'(<p) being quadratic or cubic. Here, we consider the same 
problem for the Gaussian SUSY matrix model with T = 1, where W'((p) is linear: 



W'(4>) = h{e<j>) = 60. 

We first consider the action with the external field a: 

'1 



S„ = Ntr 



-B 2 + iBe</> + if; (e ia - 1 + e) i/j 



(6.1) 



(6.2) 



Because the fermion determinant does not depend on cf) in this case, the partition function 
is essentially same as the one under the PBC (a = 0): 



e ia - 1 + e 



iV 2 



J a=0- 



(6.3) 



According to (IB.9I) in appendix [B], the Gaussian SUSY matrix model preserves the SUSY. 
Thus, the constant twist gives no effect for SUSY breaking. 



6.1 Z2-symmetric deformation 

Now the twist is deformed to (13.221) with V{e<f>) = (e0) 2 , which preserves the Z 2 symmetry 
under the flip of the signs of both and B. 

Then, the saddle point equation (I3.28P is obtained as 

dyp(y)P—= (i--)x, (6.4) 



x-y \2 c 2 

with c 2 = / dyy 2 p(y). Under the assumption that the support of p(x) is x G [—a, a], the 
solution is found as 

p(x) = — \/a 2 — x 2 , a = 2\/3, c 2 = 3. (6.5) 

67T 

Then the SUSY is spontaneously broken, since 

lim ( lim ( — tr B 2 ) ) = 1 - c 2 = -2 ^ 0. (6.6) 



a^O \ N^oo \ N 
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6.2 ^-breaking deformation 

We here consider another deformation of the twist to a linear function of (ft: a-Mr(e</>), 
corresponding to V(e<j)) = e<p, which breaks the Z 2 symmetry. (We have not considered 
this kind of deformation for the quadratic or cubic W'(4>) in the previous sections, due 
to the technical complexity.) What is interesting in this case is that, as we will see, two 
kinds of spontaneous symmetry breaking with respect to the SUSY and the Z 2 symmetry 
take place at the same time. 

The large- N saddle point equation from (13.271) 

dy P (y)P + — - - x = (6.7) 

x — y 1 — e iaCl 
with ci = J dxxp(x) is reduced to 

dyp(y)P— = i(x--Y (6.8) 



x — y 2 \ ci 

after the limit a — > 0. Assuming the support of p(x) as x G [a, b], it turns out that 
eq. (16. 8p has two solutions as 

I): a = -l, b = 3, c x = 1, (6.9) 

II): a = -3, 6 = 1, ci = -1. (6.10) 

For both solutions, p(x) takes the form 

P( x ) = ^~V( X - a )( h - x )- (6-11) 

The Z 2 symmetry is spontaneously broken at each of the saddle points. They are con- 
nected to each other by the Z 2 transformation. The semi-circle shape of p(x) is not 
influenced by changing the ^-symmetric twist to the ^-breaking one. 



6.2.1 Superselection rule 

In order to see which saddle point should be chosen, let us reconsider the saddle point 
equation (16. 7ft with keeping the linear term with respect to a. 
Then, (16. 7ft becomes 

f d yp(y)P— = i( x -L-2a) (6.12) 
J x-y 2\ c x J 

with a = jia. For the purpose of obtaining physically acceptable solutions from this, it 
is legitimate to regard the small deformation parameter a as pure imaginary rather than 
real, and a as a real quantity. Taking account of 0(a) corrections, the above solutions 
become 

I): a = — 1 + a, b = 3 + a, c\ = 1 + a, (6.13) 
II): a = s + a, b = l + &, d = -! + &. (6.14) 
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Using 



r>b /■» -y r-0 

dx dy p(x)p(y) P In \x — y\ = — dx p(x) \n(x — a) — — \ dx™ 1 

J a ^ J a 



x p{x) 



+ (^T + d ) (ci-a) + ^a 2 (6.15) 

obtained from integrating ( 16. 12ft . the planar free energy ( 13.261) can be calculated as 

F a = -ln(e 4Aci -l) 

1 1 / a + b\ 2 a 1 9 
+ 4 + 4(— j + + «(«.-)■ (6.16) 

Recall that the first term in the right-hand side of (16.161) comes from the fermion de- 
terminant (e 4aci — l)^ . So, when N is even, the first term should be replaced with 
-In \e A&Cl - l|. 

Let us consider the case N everjjj]. The expression for the free energy becomes 

I): F a = -ln|4d| + ^-2d + C(d 2 ), (6.17) 
II): F a = -ln|4d| + ^ + 2d + C(d 2 ). (6.18) 

The values of the free energy coincide up to the order O(a ). Comparing the 0(a) terms, 
we can see that I) is favored when a > 0, while II) is favored when a < 0. 
For the expectation value of B, 

lim ( lim (—ixB) J = -%c x = \~ % . !° r !L (6.19) 



a^O \ JV^oo \ N ) I +i for II) 

\ \ f as v / 

also 

lim ^lim ^trS 2 ^ ^ = -1 for I) and II), (6.20) 

implying that the SUSY is broken at either of I) or II). In this case, B plays the roles of 
two kinds of order parameters: one is for the SUSY breaking and the other is for the Z2 
symmetry breaking. 



12 For the case N odd, Z' a becomes negative at I) in (|6.13[) when a < and at II) in (|6.14|) when 
a > 0. Because Z' a is analogous to the Witten index rather than the thermal partition function, the 
negative sign of Z' a would not immediately mean instability of the system. If we consider the free energy 
associated with the absolute value of Z' a and choose the superselection sector by the value of — \&\Z' a \, 
the conclusion is same as in the case of even N discussed below. 
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6.2.2 SUSY Ward-Takahashi identity 

The SUSY breaking discussed above can also be seen as a break down of the SUSY Ward- 
Takahashi identity. For example, from the Q-transformation Q(ip<f) = —iBtfi — ipip, if the 
SUSY is not broken, the identity 

-i^ltr(S^ = ^ltr(^)\ (6.21) 

is expected to hold. 

In the presence of the external field, we compute the left- and right-hand sides: 
(l.h.s.) = -i/itr(B0 



1 

Z 



a 



d 



V (-^tr0 2 ) (e**** - 1 + ef e^K^), (6.22) 



(r.h.s.) = /-^tr(^V) 

= ±- (_i) /"^ 2 0( e -^_ 1 + e) ^-i e -ivt r (^)_ (6 23) 
Z a J 

When N is finite, Z a — > 1 as a — > 0. (See (1B.3P with #i = e > 0.) Thus, the identity 
holds as (l.h.s.) = (r.h.s.) = — 1/e. 

On the other hand, taking the large- N limit with e — * 0, and then a — > finally, we 
have 

ex(l.h.s.) -J±* M —2, (6.24, 

implying the SUSY breaking, e x (l.h.s.) jumps from —1 to —2 at JV = oo, while e x (r.h.s.) 
jumps oppositely from — 1 to at N = oo. 



6.3 General deformation 

In order to see the response of the model under further deformations, in this subsection 
we consider a general deformation of the model. The large- iV saddle point equation (13.281) 
takes the form 




(6.26) 



for x in the support of p(x), with 

eS J*VWM. (6.27) 

For two simple examples of the V(x): (one is general linear function and the other is 
cubic), large- N solutions are presented. 
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6.3.1 Example 1: general linear V(x) 
Let us consider the case 

V(x)=v + x (6.28) 

with v constant. Because (16.261) is invariant under the overall rescaling of V(x), any 
general linear V(x) can be reduced to (16. 28ft . 

Assuming the support of p(x) as x G [a, b], the two solutions are obtained: 

1 



I): a = -2 + -(-v + y/v 2 + 4), b = 2 + ~{-v + Vv 2 + 4), 



c= -(v + Vv 2 + A), (6.29) 
II): a = -2- -{v + V^TI), 6 = 2-i(*; + vV + 4), 



where 



vV + 4), (6.30) 

p(x) = ^^"^-^ (6-31) 

for both cases. These solutions reflect the spontaneous breaking of the Z2 symmetry 
again. For the expectation values of B at the saddle point I) or II), 

lim ( lim / — txB) J = -z(c - u) = -i- (-v ± Vv 2 + 4W 0, (6.32) 
q^o yv^oo \ N / J 2 V / 

1 L n2 \ \ / TV + Vv 2 + 4 



IL m „ Ufc \N tlB /J = ~ I ' 5 I • 

implies SUSY breaking. 

Since and the expectation values of B depend on v even after a — * 0, we see that, 
in the large- iV and slowly varying limit, the effect of the constant term v remains even 
after the external field turned off. 

In the sense of the critical behavior in the one- matrix model [221 [TBI El 03], from 
the semi-circle shape of p(x), the solutions represents the (2, l)-phase describing the two- 
dimensional topological gravity. 

6.3.2 Example 2: cubic V(x) 
Here, we consider 

V(x) =v + ^x 3 (6.34) 

with v constant. 

Assuming the support of p(x) as x G [a, b], the solution is given as 

a = a - ^/2a{c-a), b = a + y/2a{c-a), (6.35) 
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where a = and c are determined by 



0"(c — cr)(c — 2cr) = 2c, 

c - v = ^(9c - 5a) - ^(5c - 2a)(c - a) 2 . 



The eigenvalue distribution has the form 



P( x ) = 77— ( c - o- - x)y/Jx - a)(b - x), 

ZTTC 



(6.36) 
(6.37) 

(6.38) 



from which the model is found to belong to the topological gravity phase ((2, 1)— phase) 
without fine tuning. We will show below that the model can reach the two-dimensional 
pure gravity phase ((2, 3)-phase) by tuning v. SUSY breaking is indicated by the expec- 
tation values of B: 



lim lim ( — ti B 



lim lim ( — tr B' 

o^O V TV^oo \ N 



a — —(c — a) 

4c ; 



, \ 
1 - ca + — [c - a) 
4 ' 



(6.39) 
(6.40) 



If we denote the left-hand side of (I6.36f) as /(cr), the extrema of f(a) are 

3± V3 



6 



c. 



(6.41) 



The (2, 3)-critical point corresponds to the case that (I6.36P and (16.371) are satisfied at 
= 0* [12]. The following two critical solutions are found. 
I) = § • 3 3/4 (13 - 6^) with 



c* = 2 • 3 3/4 , a, = {VS - 1) • 3 1/4 , 

a, = -(3 - V3) • 3 1 / 4 , h = {V3 + l)- 3 1 / 4 , 

P*( x ) = 7TTr^ui( b * - x) 3/2 y/x - a*, 



2tt 2 • 3 3 / 4 



(6.42) 



II) 



3 3 / 4 (13 - 6^) with 



P*(a;) 



-2-3 3/4 , a, = -(v / 3-l)-3 1/4 , 

-(v^+l)^ 1 / 4 , 6, = (3 - \/3) • 3 1/4 , 
1 1 



2vr 2 • 3 3 / 4 



-a,) 3 / 2 v^ 



,r. 



(6.43) 



I) and II) are related by the Z 2 transformation. The fractional power 3/2 in the fall-off of 
p*(x) at the edge a* or 6* is characteristic of the (2, 3)-critical point [221 EEH1 EH EES]- The 
critical behavior describing the two-dimensional pure gravity is obtained by approaching 
6* at I) (a* at II)). 



36 



At the critical point, the expectation values of B become 



lim lim ( — txB 



lim ( lim ( — tri? 2 



^2-^3 _ 3 3 /4 



for 



I) 



2 — { II). 

10-6^ for I) and II). 



(6.44) 
(6.45) 



Similarly, higher critical points of the one-matrix models, (2, 2k — l)-critical points 
{k > 3) [221 EE31 HH US], could be accessed by the deformation with some polynomial 
V(x) of the degree 2k — 1. Our observation that tuning V(x) yields various phases of 
two-dimensional gravity manifests the fact that choosing the twist itself is a part of the 
definition of the model as stated in introduction. 



7 Summary and discussion 

In this paper, we first discussed spontaneous SUSY breaking in the (discretized) quantum 
mechanics. The twist a, playing a role of the external field, is introduced to detect the 
SUSY breaking, as well as to regularize the supersymmetric partition function (analogous 
to the Witten index) which becomes zero when the SUSY is broken. Differently from 
spontaneous breaking of ordinary (bosonic) symmetry, SUSY breaking does not require 
cooperative phenomena and can take place even in the discretized quantum mechanics 
with a finite number of discretized time steps. There is such a possibility, when the 
supersymmetric partition function vanishes. In general, some non-analytic behavior is 
necessary for spontaneous symmetry breaking to take place. For SUSY breaking in the 
finite system, it can be understood that the non-analyticity comes from the vanishing 
partition function. 

We next gave an overview of a class of large- N SUSY matrix models with the slowly 
varying parameter e introduced. The twist a is also introduced to detect SUSY breaking, 
and it is turned off at the final step. In the large-iV limit with e kept finite, there is a 
possibility that the SUSY, which is broken for finite N, becomes restored at large N. On 
the other hand, it was found that, in the large- N limit and the slowly varying limit e — > 0, 
the matrix models undergo spontaneous SUSY breaking for general superpotential. In 
the slowly varying limit, the fermion interaction terms become negligible compared with 
the term of the external field, which is the origin of the SUSY breaking. Furthermore, the 
superselection rule holds by taking the large- N limit. In the finite systems, because the 
superselection rule does not work, we saw that the expectation values of some variables 
diverges as a — > 0. However, in the matrix models, there is no such singular behavior 
thanks to the large- N limit. 

As a concrete example of the restoration of SUSY at large N, we discussed a SUSY 
matrix model with a double- well potential (the derivative of the superpotential W(<f>) 
is quadratic). In analogy with the infinite volume limit of quantum field theory, SUSY 
is expected to be broken by the instanton effect for finite N, but it is to be restored 
due to the suppression of instantons at infinite N. We obtained the large- N solutions of 
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the matrix model, and explicitly saw a number of supersymmetric behavior ensuring the 
restoration of the SUSY. 

In the large- N and slowly varying limit in the case of W'(<p) quadratic or cubic, we 
observed spontaneous SUSY breaking, and when changing the constant twist to a <fi- 
dependent function (proportional to tr0 2 ), we also found SUSY breaking. The change 
of the twist does not deform the eigenvalue distribution qualitatively, but modifies its 
precise form. Because the external field becomes dominant in the fermion bilinear terms 
in the slowly varying limit, the model will remain somewhat sensitive to the form of the 
twist even after turning it off. Thus, we can say that the model is not fully defined by the 
matrix model action alone, and the form of the twist should also be specified to determine 
the model. For the Gaussian SUSY matrix model, where W'(<j>) is linear, the constant 
twist essentially changes nothing from the model under the PBC. But, we explicitly saw 
that some ^-dependent twists cause SUSY breaking. 

For future directions, because the twists introduced in the models seem somewhat 
related to the Wilson line of an external gauge field, the models might be reformulated 
by introducing some dynamical gauge fields. Then, if the external field is dynamically 
generated by the expectation value of the gauge field, the models with various twist 
potentials V(e<p) could be unified to a single model with gauge interactions. 

It is certainly interesting to discuss whether the mechanism of SUSY breaking in this 
paper is applicable to the reduced super Yang-Mill matrix models including the IIB matrix 
model [2]. For these models, the partition function is computed in [23] to be nonzero for 
any finite N. Thus, it is hard to expect that the SUSY is spontaneously broken for finite N, 
and in particular for the IIB matrix model it would be natural to consider a possibility of 
SUSY breaking at large N, if we believe that the IIB matrix model describes our real world. 
Since such models have flat directions, for this purpose, it would be an interesting step 
to lift (some of) the flat directions slightly by introducing the slowly varying parameter 
e with keeping some of the SUSY. To carry out the analysis, the method of the Gaussian 
expansion or improved perturbation theory would be useful [2H [251 EH1 EB I2"7j . 
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A Partition function with the twist a 

In this appendix, we show that the partition function with the TBC for the fermions 
(I2.16P can be expressed in the form (12.171) . 

Let b, be annihilation and creation operators of the fermions: 

b 2 = (tf) 2 = 0, {b,P} = l, (A.l) 

and they are represented on the Fock space {|0), |1)} as 

S|0) = 0, St|0) = |l). (A.2) 

We assume that |0), |1) have the fermion numbers F = 0,1, respectively. 
The coherent states \if)), (?/>| satisfying 

&|V)=V#>, <# = <V# (A.3) 
(ip, ijj are Grassmann numbers, and anticommute with b, b*.) are explicitly constructed as 

|V>> = |0)-^|l)=e-^|0), 

<V| = (0|-(l|^ = (0|e-^. (A.4) 

Also, 

|0)= / #^#>, (0|= f 



\1) = -J (1| = J #<# (A.5) 

Thus, we can obtain 

Tr [(-e~ ia ) F e-P H ] = (0|e^ H |0) - e~ ta (l\e-P H \l) 

dipdip (e~ ia + il)i)){i)\e- pH \ip) 

= e~ ia J di)di)e eta ^{i)\e- pH \i)). (A.6) 

Since the bosonic part of H is obvious, below we focus on the fermionic part Hp = 
tfW'b. Dividing the interval /3 into M short segments of length e: (3 = Me in (1A.6[) and 
applying the relations 



(V|V) = e W 1 = J ##|^)e^(V| (A.7) 
to each segment, we have the following expression: 

Tr[(-e- ia ) F e-P HF ] = -e~ ia J njdipjd^ 



x exp 



v?=i 

M 



IVl / 

- £ Y,^ ( 

.7=1 V 



(A. 
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with 



M+l 



e ia ip 



or 



Tr [(- e ~ iQ ) F e -^] 



M 



x exp 



j^J dip j dip j 



with 



Since f1A.9h and (lA.lip correspond to (12.161) in the continuum limit e — > 0, M 
(3 = Me fixed, we find that the formula (12.1 7p holds. 



(A.9) 



(A.10) 

(A.ll) 
oo with 



B Gaussian SUSY matrix model 

In this appendix, we compute the partition function of the Gaussian SUSY matrix model 
with T = 1 under the PBC, and see that the SUSY is not spontaneously broken. From 
(13.61) and (13 .TP with a = 0, the action and the partition function are defined by 



Sq = Ntr 

zc = (-i; 



B 2 + i 9l B(f> + gi ipip 



d N 'Bd N ^d N 'ipd N 'iPe- 



(B.l) 
(B.2) 



The corresponding superpotential is given by W'(<p) = gi<p, and g\ is a (positive or 
negative) coupling constant. 

Under the normalization (13.81) . it is easy to see 



Z G = (sgnfoi)) 



N 2 



(sgn(^i)) 



N 



(B.3) 



for any N. On the other hand, the partition function at large N is evaluated by a saddle 
point configuration of the eigenvalue distribution of 0: 



PG{X) 



N 



trS(x — 0) 



(B.4) 



as 



+ O{N 



dxdyp G (x)pG(y)Pln\x -y\- / dx-x 2 p G (x) 



(B.5) 
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Cn is a numerical factor whose large- N behavior is fixed below. 

Using the standard method in [19J, it turns out that Pg(x) has a semi-circle shape in 
the support x G [—2,2]: 



PG[X) 



2tt 



X/ 



The exponential factor of ( 1B.5j) can be explicitly computed to give 



Z G = (sgn(5f 1 )) 7V C N exp 



~N 2 + O(N ) 



(B.6) 



(B.7) 



Comparing this with (1B.3I) . the large- N behavior of Cjv is determined as 

"3 



Cn = exp 



-N 2 + O{N ) 



(B.8) 



Since the ^-integrals in (IB .20 yield the delta-function of B, we can see that all the 
expectation values of -^tr B n (n = 1, 2, • • • ) vanish 



N 



-MB n 



(B.9) 



which implies that the SUSY is not spontaneously broken. 



C Proof of (14^1 ) 



In this appendix, we give a proof of ( 14. 9p . which serves as a strong evidence for the 
restoration of SUSY in the large- N matrix model with a double-well potential discussed 
in section 14.1.11 

First, (14.91) holds for n — 1,2 as seen by the explicit computation. 
Next, let us fix an arbitrary n > 2. In order to prove (14. 9 j) by induction, it is sufficient 
to show (^tr B n+1 ^ = under the assumption 

itrS*\=0 (k = l,2,-.. ,n). (C.l) 

In this appendix, for notational simplicity, we write neither the symbols lim a ^ > h m 7v-»oo 
nor the suffix a of the expectation values. 

For our purpose it is useful to notice a Schwinger-Dyson equation for k > 1, m > 0, 

1=0 \ ' ^ 

- /^tr {B k+1 W'^T) \ - ,7itr !U'(0) m+1 ) ) , (C.2) 
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derived by taking the large- N limit of the identity 

N 2 



= f d N2 (j>d N2 Bd N ^d N ^y^ti [BH^W'm e 



(C.3) 



where t@ is an appropriate basis of the N x N hermitian matrices. Then, by the assump- 
tion, eq. (IC.2I) is reduced to 



1 

N' 



trB 



n+l 



-/< -tx{B n W'{4>)) 



(C.4) 



for k = n, m = 0, and 
1 



N 



tr (B k+1 W'{<f)) m ) 



/( l-M{B k W\<P) m+1 )) + (^tT{B k ~ 1 W'{<P) m )) (C.5) 



N 



for 1 < k < n — 1, m > 1. 

By using these equations repeatedly, we can lower the power of B in (^tr B n+1 }, and 
eventually (-^trl? n+1 ) is expressed as a linear combination of \-htT W (0) m ) for m = 
0, • • • , n + 1. The coefficients of the linear combination are given by the combinatorics 



of how many times we use the first and the second term on the right-hand side in ( 1C.5I) . 
Here, it is convenient to associate a monomial s k+1 t m with (^tr (B k+1 W'((f)) m ^y Then, 
( 10.4ft and (10. 5ft tell that we can make the following replacement: 



,n+l 



s k+i t m 



s n+1 x (-i)-, 



s k+l t m x | , { 



t 1 



(C.6) 



Then, (^tr B n+1 } is eventually associated with s°-terms by making these replacements 
repeatedly starting from s n+1 . It amounts to the s°-terms of 



— is n t x 



m=0 x ' 



-is n t ■ 



1 + V 

s S A 



(C.7) 



which can be expressed as 



s n t 



c 



ds 1 



2m s 1 + it - \ 



ds 



C 



2ni s 2 + ist — 1 ' 



{CI 



where C is a sufficiently large circle around the origin. It is easy to compute this integral 

as 

-it 



-it + V4^¥ 



n+l 



-it - V4^t2 



n+l 



(C.9) 



(It turns out to be a polynomial of t of the degree n + l.) 
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Notice that t corresponds to W'{4>). Using 



-tr w'(<py 



dx p(x)(x 2 + jj, 



dt 



2\m 



2 2tt 



we find 



— tr B n+l 
N 



2 dt 

— t 
-2 2tt 



-it + a/4 - & 



n+1 



-it - V4 - t< 



n+1' 



(CIO) 



(c.n; 



It is straightforward to see this integral vanishes for all n > 2. Thus, the induction 
argument completes to prove (14. 9p . 
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